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Q , We present results of 1.3 mm interferometric and single-dish observations 

of the center of the L1641-N cluster in Orion. Single-dish wide- field contin- 
'■ uum and C0(2-l) observations reveal the presence of several molecular outflows 

^ ' driven by deeply embedded protostellar sources. At higher angular resolution, 



O 



the dominant millimeter source in the cluster center is resolved into a pair of 
protostars (L1641-N-MM1 and MM3), each driving a coUimated outflow, and 
a more extended, clumpy core. Low-velocity CO line-wing emission is widely 
! spread over much of the cluster area. We detect and map the distribution of 

several other molecular transitions (^^CO, C^^O, ^^CS, SO, CH3OH, CH3CN, 
^ i and OCS). CH3CN and OCS may indicate the presence of a hot corino around 

^ ; L1641-N-MM1. We tentatively identify a velocity gradient over L1641-N-MM1 in 

CH3CN and OCS, oriented roughly perpendicular to the outflow direction, per- 
haps indicative of a circumstellar disk. An analysis of the energy and momentum 
load of the CO outflows, along with the notion that apparently a large volume 
fraction is affected by the multiple outflow activity, suggests that outflows from 
a population of low-mass stars might have a significant impact on clustered (and 
potentially high-mass) star formation. 
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mation — turbulence 
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Introduction 



Most stars are born in clustered environments (e.g.. lLada &: Ladall2003l ) and their for- 
mation may be significantly different than in isolated regions such as the well studied Taurus 
cloud. For example, the shape and extent of prestellar cores and their subsequent evolution 
may be affected through close encounters with other cluster members or the effect of ion- 
izing radiation from massive cluster stars. Protostellar outflows, known to accompany star 
formation, might also play a different role in clusters. Whereas in the isolated case much 
of the outflow energy and momentum will likely just be dumped somewhere else in or even 
outside the cloud, multiple, randomly oriented outflows in a cluster-forming clump might 
impact a larger volume, significantly affecting or even hindering further star formation. 

We present new data on the outflow activity in the L1641-N region in the Orion A 
giant molecular cloud. As a moderately sized, relatively nearbj0 cluster, L1641-N allows for 
detailed observations, in a relatively simple example, of the complex processes involved in 
cluster formation. 

The first detection of b i polar CO outflow activity in L1641-N, ce ntered on IRAS 053 38—0624, 
was found bv iFukui et al.l (Il986l ) and subsequently followed up by iFukui et al.l (llQSSi ) and 
Wilking et al.l (jl990l ). Near infrare d imaging revealed the presence of an embedded small 
group of about 20-25 yo u ng st ars (IStrom et al.l Il989al : I Chen et al.l Il993l : iHodapp fc Deane 
19931 ). iHodapp fc Deand (119931 ) estimated that the cluster grew at a rate of 2-3 stars every 
10^ yr over about 10^ yr. 

The CO outflow maps have generally been interpreted as evidence for the presence of 
one, roughly north-south oriented outflow, driven by one protostellar source seen as the mid- 
to far-infrared source IRAS 05338—0624. Ground based M-band observations have revealed 



a source Nl seen at 5 /zm (and longward) only (IChen et al.l 119931 ): thi s source is also seen 



as a compact raillimeter continuum so urce in interferometric images (IWilking et al.l Il989 



McMuUin et al.ll 19941 : IChen et al.lll995l). but app e ars not to be the dominant source at least 



at mid-infrared wavelengths ( Stanke et al. 1998 : Ali fc Noriega- Crespol 



2004, AN04 in the 



following). Moreover, interferometric CS and SiO observations (iMcMuUin et al.lll994j ) sug- 
gest a south-west to north-east oriented outflow at the very cluster center, as opposed to the 
predominantly north-south orientation seen on larger scales in CO. IStanke et al.l (|l998l ) have 
proposed that there are indeed 2 independent outflows: a south-west to north-east oriented 
flow, and a large-scale, north-south outflow, driven by the brightest 10 /xm source in the field. 
The south-west to north-east flow appears to be driven by the compact millimeter continuum 



^we assume the canonical Orion distance of 450 pc 
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source. Further optical and near-infrared narrow-band imaging searches for Herbig-Haro ob- 
jects and shocked H2 emission have provided evidence that the L1641 - N cluster indeed is 



Mader et al. 1999 



the center of mult i ple outflows (IDavis et al.l Il995l : iReipurth et al.l Il998l : IStanke et al.l Il998 



Stanke et al.l 12001 l2002l . (hereafter SMZ2002)), some of which extend 



over several parsecs. 

In this contribution we present new 1.3 mm wavelength single-dish and interferometric 
observations of the L1641-N cluster area in order to shed further light on the protostellar 
outflow activity. We will first describe the observations (Section [2]), identify outflows and 
their driving sources in Section [3l and provide a summary and conclusions in Section |H 



2. Observations & data reduction 
2.1. IRAM 30 m 

We obtained 1.3 mm dust continuum maps in the course of a wider survey over several 
observing runs between 1999 and 2002, using the 37 channel MAMBO bolometer array at the 
IRAM 30m telescope. The data reduction was done with MOPsfl and followed standard 
bolometer reduction procedures, including baseline subtractions, despiking, correction for 
athmospheric extinction using the results of skydip measurements, and flux calibration using 
planet maps. Correlated sky-brightness variations (sky-noise) were removed in an iterative 
process, using the resulting mosaic of the respective previous iteration as input source model 
to the actual sky-noise reduction. Thereafter, the single beam intensity distribution was 
restored from the chopped dual-beam data and subsequently combined into the large mosaic. 



The 9 element heterodyne receiver array HERA (jSchuster et al.l |2004J ) was used in 
November 2003 to obtain C0(2-l) maps of the L1641-N region. The cluster center and 
the area around the driving source of Outflow 51 in SMZ2002 were each covered with a fully 
sampled, 66"x66" raster map. In Addition, a ~10'x9.'5 area was mapped in the coarsely 
sampled spectral line On-The-Fly mod^. Data reduction included low-order baseline sub- 
traction and combining the spectra (raster and OTP maps) into one map. 

To convert CO brightness temperatures to molecular gas column densities, we assumed 
the CO to be optically thin, in LTE at a temperature of 30 K, and a CO abundance of 
10~^ of the H2 abundance. We derived gas masses as a function of velocity for each channel. 



^developed and maintained by Robert Zylka, IRAM Grenoble, France 

■^sampling at 7.6"; see HERA manual for details: http: / /www. iram.es /IRAMES / mainWiki/HeraWebPage 



-4- 



Kinetic energies and momenta were derived as i?kin = ' ~ ^ccn)^ and P = 

M{v) ■ {v — Veen), with t>cen = 6.8 km/s. A characteristic timescale was derived by dividing 
the maximum distance / at which high-velocity emission was seen by the maximum CO 
velocity observed. Results for the various outflow lobes identified from the data are listed in 
Tab.E 

We do neither attempt to correct for optical depth effects in calculating the CO outflow 
masses nor for outflow gas at radial velocities equal to the ambient cloud velocities; similarly, 
we do not correct the velocity of the gas for projection effects. Outflow masses, kinetic 
energies, and momenta are therefore strict lower limits. 

2.2. SMA 

We used the Submillimeter Array (SMA) on Mauna Kea, Hawai'i, to obtain high res- 
olution 230 GHz images of the center of the L1641-N cluster. Two tracks were taken, one 
in the compact (Dec. 19, 2004, with antennas 1-6; antenna 8 was flagged) and one in the 
extended configuration (Dec. 27, 2004, with antennas 1-6 and 8). The two configurations 
yielded projected baselines ranging from about 13 to 220 m. Double sideband receivers were 
used, with the CO (2-1) line centered in the upper sideband (correlator chunk 14). 230 GHz 
zenith opacities were between 0.1 and 0.06 during the first track, with fairly stable phases; 
system temperatures were between 200 and 250 K throughout most of the track. During the 
second track, the opacity was initially high, ~0.2, but dropped to ~0.08 later on; similarly, 
phase stability improved during the track. System temperatures were around 150 K. 

Data calibration and editing were done within the IDL MIR package. The phase center 
for L1641-N was set to a=5^36""18.'^8, 5=-6°22'10"(J2000). Observations of Uranus and 
Callisto were used for bandpass and flux calibration. Flux calibration was obtained using 
model predictions for Callisto's flux depending on baseline length as provided on the SMA 
webpage^. Observations of J0423— 013 and J0607— 157 were taken every 20 min to allow 
for complex phase calibration. Continuum emission was reconstructed from line-free parts 
of the 2x2 GHz range covered by the SMA receivers. 

Maps were constructed and deconvolved within MIRIAD. For the emission line maps, 
continuum emission was first subtracted from the emission line visibilities. The continuum 
and ^^CO line SMA data were complemented with short spacing visibilities obtained from 
the IRAM 30 m mapping described above. Maps were created using a natural weighting 



^http:/ /smal.sma. hawaii.edu/ 
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scheme, and CLEANing was performed, stopping after new CLEAN components were at 
about 60% of the rms noise level. The synthesized beam measures l'.'7xl'.'2 at a position 
angle of 75° east of north. 



2.3. Astrometry 

In this paper we use a number of different data sets, taken over a wide wavelength 
range: near-infrared (i^'-band) data from SMZ2002, ISO data published by AN04, and new 
single-dish and interferometric millimeter wavelength data. Good relative astrometry is 
crucial in understanding what happens in a clustered environment such as L1641-N. The 
various datasets used in this paper were registered with respect to each other as follows: the 
infrared K-hand data from SMZ2002 were recalibrated over the field of interest using 2MASS 
coordinates for stellar sources. The ISO maps were registered to the near-infrared using a 
number of sources seen in the near- and mid- infrared maps. The 1.2 mm dust continuum 
map was registered with respect to the near-infrared map using one star seen in both maps, 
located 1'.5 south-east of the cluster center. The interferometer map was not corrected, 
because it is phase referenced to quasars, yielding high positional accuracy. 



3. Young stellar objects and outfiovi^s in L1641-N 

Figure [1] shows an overview of the field around the L1641-N cluster in 1.3 mm dust 
co ntinuum (contours) a nd 2.12 fim narrow band (greyscale: H2 f =1-0 S(l) line + continuum; 



see 



Stanke et al.l JlQQsh . SMZ2002). Fig.[2]shows the CO J =2-1 line- wing emission (blue/red 



contours) . 



3.1. Wide- field mapping 



The wide- field 1.3 mm dust continuum map taken with MAMBO shows a number of 
compact features in the center of the field and more diffuse emission towards the south-east 
(FigJll). This i s remi niscint of the head-tail s hape of the C^^O cloud core as mapped by 



Reipurth et al.l (Il998l ) (their Fig. 17a) (see also iTatematsu et al.lll993l ). 



The CO line-wing emission maps (Figs. [3] and E]) reveal the presence of a number of 
well-defined outflow lobes. The dominant feature at red-shifted velocities is a ~lf5 broad, 
collimated lobe extending over 7' south of the L1641-N cluster center to the edge of the area 
mapped in CO (labeled R-S in Fig. [2]). It is also seen as a series of H2 shocks reaching even 
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further south (jStanke et aLll2000l ). Together with a chain of optical Herbi g-Haro objects 
north of the cluster (HH 306-HH 310: iReipurth et alJEoOsI : iMader et al.lll999h these features 
form a bipolar giant outflow stretching over several parsecs on both sides of the cluster. 
Blue-shifted CO emission is also seen north of the cluster center in a 'Y'-shaped feature. 
The channel maps in Fig. [3] show that at high velocities there is a collimated lobe extending 
NE of the cluster center (B-NE), which is joined by an apparently also collimated N-S feature 
(B-N) at intermediate velocities. At low velocities, these features blend with a broad patch 
of CO emission north of the cluster center. Broadly distributed red-shifted emission is also 
seen to the south-west and west of the cluster center (Fig. [5]). High- velocity emission seems 
to be restricted to a small lobe close to the cluster center, forming a counter-lobe to B- 
NE. The MAMBO dust continuum map shows a bright, compact but resolved peak at the 
center of these CO outflow lobes, on top of a north- east to south-west oriented ridge, which 
corresponds to the HCN and HCO"*" ridge found by lFukui et al.l (Il988l ). 

Another well defined, highly collimated outflow lobe (R-SW) is seen to the west of the 
large N-S lobe R-S at a position angle of about 28° east of north. The CO channel maps 
(Fig. ED reveal a distinct velocity structure in this lobe, with the highest velocity emission 
being found at the south-western tip, and lower velocity emission extending ever further to 
the north-east. In fact, close to the ambient CO velocity (Fig. Hj) this lobe structure can 
be clearly followed up to where it intersects with the large, redshifted N-S giant flow lobe. 
The MAMBO continuum map also shows a faint ridge alo ng the northern p ortion of this 
CO lobe. We suggest this is outflow entrained dust (c.f., IChini et al.l I2OOII ). Converting 
the CO flux densities to broad-band flux densities yields flux densities of the order of 2- 
5 mJy/beam (using a conversion factor of 8.6 Jy/K as given on the IRAM webpages and a 
MAMBO bandwidth of 100 GHz estimated from the transmission plots shown on the MPIfR 
MAMBO webpages), much less than the 30-50 mJy/beam seen in the MAMBO map in the 
ridge. The R-SW lobe does not have any obvious blue-shifted counterlobe. No H2 shocks 
have been found that could be associated with this lobe by SMZ2002. 

A bipolar CO outflow is centered on a=5*^36™24.^8, 5=-6°22'42". It corresponds to 
[SMZ2002] Outflow 51. The outflow is centered on a north-west to south-east oriented 
millimeter continuum ridge and an ISO mid- infrared source (source 10 of AN04). Its 
blueshifted lobe is seen over a large distance from the so urce as a sequence of large optical 



HH-objects (HH 301 an d HH 302) and infrared shocks (IReipurth et al.l 1 19981 : iMader et al. 



I999I : IStanke et al.ll2002l ). No such features have been found at comparable distances from 



the source on the redshifted side. We speculate that the flow might have a continuation on 
the redshifted side in the outflow lobe R-SW, for which we cannot see any other plausible 
driving source. This would imply a significant change in outflow direction, however, and 
cannot be due to precession as only the redshifted lobe shows this bend. We noted that at 
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velocities close to the CO ambient cloud velocity the R-SW lobe can be traced back to where 
it intersects with the R-S lobe. This is also about where the R-E lobe of [SMZ2002] Out- 
flow 51 would intersect with R-S. We therefore speculate that Outflow 51 interacts with 
R-S, its outflowing gas being accelerated in a southerly direction, thus creating the bend 
in the redshifted flow lobe . Deflection of a jet with a side wind is a candidate mechanism 



to explain bent jets (e.g., iFendt &: Zinneckei 



IQOsh. particu larly in the case of 'C'-shaped 



symmetry about the driving source (e.g., iBally et al.l l2006l ). Observations of jets of this 
type, as well as theoretical and experimental studies, demonstrate that jets can survive (and 
remain coUin aated) while interacting with relatively tenou s winds such as expanding HII 



regions (e.g., iMasciadri fc Ragall200ll : iLebedev et al.ll2004l ) over a large portion of the jet 
length. In L1641-N in contrast, the R-SW lobe seems to interact with R-S only over a short 
total length, and it seems that R-S is heavier than R-SW, as R-S does not show any signs of 
bending resulting from an interaction with R-SW. That a lighter jet might be able to survive 
a collision with a heavier flow remains to be shown, but it seems that a collision of two jets 



and redirection does not necessarily imply a disruption of the flow(s) (ICunningham et al. 



20061). 



In additio n, we identify a s mall, bipolar CO outflow centered on a near- to mid-IR 
(source N31 of IChen et al.l Il993l . ISO source 13 of AN04) and 1.3 mm continuum source 
at a=5^^36'^22.^0, (5=-6°23'28". The blueshifted lobe is associated with two faint H2 shocks 
which were not identified by SMZ2002. 



Finally, the Strom 11 group of embedded stars (IStrom et al.lll989bl : IChen et al.lll993l ) 
is seen superimposed on a roughly circular millimeter continuum clump of ~1' diameter. A 
faint, narrow lobe of blueshifted CO emission (B-SE2, seen best in the channel maps in Fig. 
[3l close to the ambient velocity) might have its origin in the Strom 11 group; some more 
faint, blueshifted emission is seen in Strom 11 and to its south. 



3.2. The central area: SMA observations 

3.2.1. Continuum 

Figure [6] shows the 1.3 mm dust continuum at high angular resolution (SMA and 
MAMBO data combined). We detect a bright, compact source centered at a =5^ 36™18.'^780, 



(5^— 6°22'10'.'37 (J2000), in very good agreement with the position found by IChen et al. 



(119951 ). We will refer to this source as L1641-N MMl in the following. The visibilities de- 
crease from 0.3-0.4 Jy at short spacings to 0.1-0.2 Jy at large uv-distances, so the source is 
resolved. A gaussian fit to the radially averaged real continuum visibilities implies a FWHM 
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size of I'/O ± 0.1. From elliptical gaussian fits in the image plane we find a size of 2'.'lxl'.'5 
(P.A. 67°), corresponding to a deconvolved size l'.'2x0'.'8 (P. A. 56°). Fitting the combined 
MAMBO+SMA data at fiuxes >0.02 Jy, gives a size 2'.'5xl'.'7 (P.A. 67°), corresponding to 
l'.'9xl'.'3 (P.A. 62°) deconvolved. The image shows a tail of emission running at a similar po- 
sition angle away from the source, so the P.A. derived from the gaussian fit might be biased 
towards this direction. As the deconvolved source size is significantly smaller than the beam 
size, the result is likely subject to large errors. Still it seems that the continuum emission is 
not oriented perpendicular to the outfiow direction inferred from the CO observations (see 
below). As the emission of MMl furthermore blends into more extended core emission, we 
suggest that MMl traces an envelope structure rather than a disk. The total integrated 
fiux of 0.54 Jy (as measured from the gaussian fit) corresponds to a circumstellar gas mass 
of about 1.6 Mf7), assuming a dust temperature of 20 K and a dust opacity coefficient of 



0.01 cmV^ flMotte et al.lll998f ). 



To the south of MMl a number of extended features form a clumpy core. The most 
prominent of these are two clumps to the south of L1641-N MMl (L1641-N MM2 and 
L1641-N MM4 in the following). There appears to be an extension of L1641-N MMl in 
a SW direction, forming a bridge of emission between L1641-N MMl and L1641-N MM2. 
Another extension to L1641-N MMl is seen to protrude towards the ESE, which appears as 
a separate source in some of the observed molecules (see below) ; we will refer to this feature 
as L1641-N MM3. 

Protruding to the NE of L1641-N MM3 there appears to be a clumpy ridge of emission 
parallel to the small H2 jet (feature C in Fig. [6]). This ridge might outline the wall of a 
cavity produced by the jet. Finally, some more clumpy emission is seen to the south and 
west. Generally, there is a good anticorrelation between the millimeter continuum emission 
and the diffuse infrared emission, indicating that the dust is in the foreground of the diffuse 
emission. 

Notably, the single-dish continuum peak does not coincide with the position of MMl; 
although pointing errors might affect the single-dish peak position, the offset seems to be 
too large; note also that the astrometry of the single-dish map has been corrected using a 
near-infrared source, the position of which should be accurate to better than 1 arcsecond. 
Thus the single-dish continuum peak seems to trace the clumpy core south of MMl rather 
than MMl itself. 
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3.2.2. CO line-wing emission 

Figure [7] shows the spatial distribution of CO hue- wing emission over the field observed 
with the SMA. The emission has been split up into a low- and high-velocity component 
(left /right panel, respectively). At high velocities, there is evidence for the presence of a well 
coUimated (though faint), NE-SW oriented molecular jet. Apparently, it is not centered on 
the bright MMl continuum source, but on MM3. The red-shifted lobe of this collimated fiow 
is more pronounced and seen over a larger distance from the driving source. It appears to 
terminate in a shock feature seen in various other molecu lar lines (SO and CH3OH, see below; 



CS and SiO fIChen et all fll996[ ): iMcMuUin et all t994 \). The blue-shifted lobe appears to 



be more spurious. Some blueshifted high-velocity emission seems to be present also towards 
the prominent H2 bow shock Al seen in the NE corner of Fig. [71 

At lower velocities, CO emission is spread over a wide area. Blue-shifted emission is seen 
virtually over the entire north-eastern quadrant, while red-shifted CO prevails towards the 
south and west of the cluster center. Centered on MMl there appears to be a narrow, hollow 
cavity-like structure, particularly clear in the blueshifted lobe. The NE continuum ridge is 
found just south of the southern edge of the blueshifted CO cavity, apparently outlining the 
wall of the outfiow cavity. The small H2 jet formed by the knots labeled as C in Fig. [6] 
appears to run right through the CO cavity; however, it cannot be excluded that this is due 
to residual error in the astrometry for the H2 image. 

The data suggest that there are two almost parallel (in projection) outfiows in the cluster 
center. One is marked by the bipolar, high-velocity CO jet and apparently driven by MM3. 
The other one is indicated by the presence of the bipolar cavity centered on MMl. 

It is straightforward to relate the single-dish B-NE lobe and its short, red-shifted coun- 
terlobe to the fiows seen at high-resolution with the SMA. The H2 bow-shock Al coincides 
with the point out to which high-velocity blue-shifted CO is seen in the single-dish data; 
however, at lower velocities, this collimated outfiow lobe extends further out to the north- 
east. Remarkably, the axis of none of the two fiows seen at high resolution with the SMA 
really points in the direction of the Al II2 bow shock. We suggest that it is MMl, which 
creates the Al bow shock and the bulk of the single-dish B-NE lobe, simply because it seems 
that MMl is the more massive object. 

The molecular outfiow from L1641-N has generally been interpreted as a single, wide- 
angle fiow. Our SMA and IRAM 30 m CO observations show that this is not true. There 
is clearly outfiow along a south-west to north-east axis, marked by the fiows seen at high- 
resolution with the SMA and the B-NE lobe. In contrast, the R-S CO lobe as well as near- 
infrared and optical shocks indicate outfiow along a N-S axis. The coexistence of outfiow 
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features along both axes within the central arcminute around the cluster center (CO lobes, 
optical and infrared H2 shocks) rule out the possibility that an initially N-S oriented flow 
could have changed its orientation to a NE-SW orientation, e.g., as a result of an interaction 
of its driving source in a close encounter with another object. Close inspection of the 30 m 
CO channel maps also indicates that the N-S outflow axis, defined by B-N and R-S, is shifted 
by several arcseconds to the east of MM1/MM3. We therefore conclude that there are indeed 
3 CO outflows from the cluster center; two, driven by MMl and MM3, seen with the SMA 
along a NE-SW or ientation, and a N-S oriented giant outflow. As suggested already by 



Stanke et al.l (119981 ) the dominant mid-infrared source in the cluster as seen in ground-based 



and ISO observations (AN04 source 18), located about 9" to the east of MM1/MM3, appears 
to drive this outflow. 



3.2.3. Other Unes 

The SMA provides a frequency coverage of 2 GHz in each of the upper and lower 
sidebands and allows several lines to be observed simultaneosly. Centering the C0(2-l) 
line in the USB yields the isotopomers ^'^CO and C^^O in the LSB. In addition to the CO 
lines, we find some CH3OH transitions, ^^CS, OCS, SO, and the first few transitions of the 
CH3CN (12-11) series. Fig. [S] shows a compilation of spectral line maps, with emission close 
to the rest velocity (5-8.6 km s~^) and line-wing emission (8.6 to 13.4 km s~^: red; 0.2 to 
5 km s~^: blue) plotted separately, and Tab. [2] gives a list of detected lines. Besides the 
lines shown in Fig. [8] there is a marginal detection of the HNCO 10(0,10)-9(0,9) transition 
at 219.798 GHz, which however is too faint for mapping and will not be discussed further. 

^^CO and C^^O trace a more extended structure (the former clearly affected by miss- 
ing short spacing information), with some emission associated with MMl and MM3. The 
brightest peak in ^^CO and C^^O is significantly offset towards the SW of MMl. MM3 is 
coincident with the second brightest ^^CO peak and a faint C^^O peak; the second brightest 
C^^O peak is seen towards the NE continuum emission ridge. Another C^^O peak i s foun d 



some 17" north of MMl, roughly coinciding with the CS peak 3 found bvlChen et al.l (119961 ). 



Besides the CS peak 1, none of the other CS peaks found by I Chen et al.l (119961 ) is seen in 
^^CO or C^^O. ^^CO shows clear signs for emission in the line-wing bins, in a bipolar fashion 
around MMl; particularly the blue-shifted emission corresponds well to the structure seen 
at the base of the cavity seen in ^^CO at moderate velocities. No ^^CO line- wing emission is 
obviously related to MM3, in agreement with our suggestion that MMl is the more massive 
outflow. 



The other molecules show pronounced differences in their distribution. SO emission is 
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found towards MMl, MM3, and an extended feature to the SW of the cluster center. Besides 
two compact features corresponding to MMl and MM3 there is also some more extended 
emission surrounding these sources, which is seen in the central velocity bin as well as in 
the line- wing bins, mostly to the east of MMl. In the integrated spectrum the blue-shifted 
wing of the SO line is much more pronounced than the red wing, whereas in the map, the 
red wing seems to dominate, particularly due to strong emission seen towards the SW shock; 
some more red-shifted SO emission is seen between the cluster center and the SW shock. 
The SO channel maps (not shown here) suggest that most of the blue-shifted SO emission is 
in extended, low surface brightness emission towards the blue-shifted lobe of the flow from 
MMl, which is too faint to show up in the map in Fig. [HI 

CH3OH emission is seen mainly towards MMl, and towards MM3 in the 8(— 1,8)- 
7(0,7) E line (and maybe marginally in the 8(0,8)-7(l,6) E line). Around MMl and MM3 
the emission in the 8(— 1,8)-7(0,7) E line appears to be somewhat extended. The SW shock is 
seen in the 8(— 1,8)-7(0,7) E and 8(0,8)-7(l,6) E lines, with a clear redshifted component in 
the 8(— 1, 8)-7(0,7) E line: it exte nds out to veloci t ies ar ound 15-16 km/s, consistent with the 



results of iMcMulhn et al.l (119941 ) and lChen et al.l (119961 ) (CS peak 2). The CH3OH 15(4,11)- 



16(3,13) line is seen as a faint, compact feature towards MMl only; the 3(— 2,2)-4(— 1,4) line 
in addition is marginally seen in MM3. 

^^CS appears to be spread out over a similar area as C^^O, but is only seen very faintly 
in the map. OCS is only seen as a compact feature coincident with MMl. CH3CN is clearly 
detected as a compact feature coincident with MMl and marginally towards MM3 (Fig. [HI 
top). 

The SO, CH3OH 8(-l,8)-7(0,7), CH3OH 8(0,8)-7(l,6), and ^^CO line show significant 
line wing emission in spectra integrated over the area covered by the SMA primary beam; 
faint wing emission might also be seen in C^^O. The wing emission is more pronounced on the 
blue-shifted side, but fainter red-shi fted emission is also pr esent. Single dish C^^O, SO and 



CH3OH spectra (CSO, ~30" beam, iMcMuUin et al.l (119941 )) tend to show more prominent 



emission in the red wing, suggesting that there is a significant extended component which 
is resolved out by the interferometer. The remaining lines (CH3OH, CH3CN, ^'^CS, OCS) 
as well as the cores of the lines associated with wing emission are resolved, with FWHMs 
around 3 km/s. 

MM2 and MM4 are not seen in any of the lines. 

Complex organic molecules are thought to form in the dense, dusty environments around 
young protostars as the pro ducts of grain surface chemistry evaporate and react with each 



other (iBottinelli et al.l 120071 ). Whereas CH3OH is clearly associated also with the outflows. 
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this seems not to be the case for CH3CN and OCS in L1641-N. We have used the XCLASS 
extension to the CLASS software (P. Schilke, pers. comm.) to simulate the CH3CN spectrum 
to estimate the temperature, column density, and size of the emitting region (Fig. [HI bottom; 
we used a Jy to K conversion factor of 12, corresponding to an observed convolved source size 
of 1'.'4, i.e., unresolved by our beam; this is a vahd assumption, as the line fluxes integrated 
over the source are found to correspond well to the peak fluxes per beam). The detection 
of the K = 4 and 5 lines {E^ = 183 and 247 K) suggest a temperature of more than 100 K, 
and their brightness with respect to the lower K lines limit the temperature to less than 
200 K. The K = 0, 1,2, and 3 lines are of similar brightness, which suggests that they are 
somewhat optically thick (for the adopted temperature range, the K = 3 line {E^^ = 133 K) 
should be clearly brighter than the K = 2 line if optically thin, due to its high statistical 
weight). The observed line intensities then imply a CH3CN column density of the order 
of a few times lO^^cm"^ and a source size of the order of 0'.'4, i.e., 150-200 AU. This is 
significantly smaller and somewhat cooler than CH.sCN emi s sion r egions found in hot cores 
in massive star forming regions (e.g., W3(H20), IChen et al.l (120061 )) but is is comparable to 
the hot 'corinos' which have now been found in a handful of other low-mass protostars (e.g.. 



BottineUi et al.ll2004l : |j0rgensen et all 120051 ). 



3.2.4- Evidence for a rotating disk in MMl 

In Fig. [To] (left panel) the positions of the emission maxima, as derived from gaussian 
fits, are plotted for individual 1.2 km s~^ wide channels for the CH3CN K=0, 1, 2, and 3 
lines and the OCS line. Only channels in which emission above 5(T is seen are used. We 
chose to use only CH3CN and OCS lines in this plot as it is these molecules where we see 
the least contamination by outflow related emission. 

Overall, the positions of the emission maxima define a south-east to north-west oriented, 
elongated structure (position angle ~ 120°). This is unlikely to be the result of positional 
measurement errors, as this would tend to produce a distribution which is elongated along 
the beams major axis, which is at a position angle of about 75°. There is a clear trend for 
the lower velocity center positions to cluster to the south-east of the MMl central position, 
and the higher velocity positions to be found to its north-west. To illustrate this further, 
the right panel of Fig. [TU] shows the channel velocities as a function of the offset along the 
disk plane. Although the positional errors are larger than the offsets, a trend for velocities 
increasing from sout-east to north-west is consistently seen in all transitions. Clearly, better 
data are needed to confirm this observation. 

We tentatively suggest that the positional arrangement of the emission maxima in- 
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dicates a flattened structur e around the protostellar s ource embedded in MMl; a similar 
result has been obtained by ICesaroni et al.l (119971 . Il999l ) for the intermediate- mass protostar 
IRAS 20126+4104, and has been interpreted as evidence for a circumstellar disk. The ve- 
locity gradient hints at rotation, with a rotation axis being roughly parallel to the axis of 
the outflow from MMl. This would imply that the CH3CN and OCS emission arises from 
grain mantle evaporation in a heated disk rather than the dense, infalling envelope. How- 
ever, higher angular resolution observations resolving the source is required to discriminate 
between these scenarios (is it spherical or flattened?) definitively. Also, neither the angular 
resolution nor the velocity resolution nor the sensitivity of our data are good enough to 
derive estimates of the central object's mass or the size of the disk. 



4. Summary and Discussion 

We have presented new 1.3 mm dust continuum and molecular line observations of the 
L1641-N cluster area. Single-dish wide field CO mapping revealed the presence of several 
blue- and red-shifted outflow lobes, which can be grouped into at least 4 bipolar flows: a 
north-south oriented giant outflow B-N/R-S, driven most likely by the AN04 ISO source 
18; a shorter north-east to south-west oriented outflow B-NE/R-W driven by the dominant 
millimeter continuum source in the cluster center; another parsec scale outflow B-E/R-E/R- 
SW corresponding to the H2 flow 51 of SMZ2002, driven by an embedded ISO and millimeter 
source east of the cluster center; and a small flow B-SE/R-SE driven by a near-infrared star 
south-east of the cluster center. Moreover, some CO outflow activity might also be associated 
with the Strom 11 group further south-east of the L1641-N cluster center (B-SE2). 

The SMA observations show that there are in fact two deeply embedded sources in the 
very cluster center, MMl and MM3, each of which drives a CO outflow, which together 
make up the B-NE/R-W outflow: MMl is at the center of two cavity-shaped outflow lobes 
seen at moderate CO velocities, and MM3 drives a collimated, faint jet seen only at high 
velocities. The flow from MMl carries more mass as only this flow is seen in the optically 
thin ^'^CO. We find ^^CO, C^^O, and ^^CS emission tracing the core from which MMl and 
MM3 are forming, CH3OH and SO emission tracing the protostellar sources MMl and MM3 
as well as their outflows, and OCS and CH3CN emission from MMl. Analysis of the CH3CN 
J = 12 — 11 {K = 0...5) lines suggests that the emitting gas traces a compact (100-200 AU), 
warm (100-200 K) region deep inside the MMl envelope. This could be the warm, innermost 
portion of the infalling envelope, i.e., a hot corino, as seen in other low- mass protostars. 
Alternatively, it may be the the heated surface layer of a disk, as suggested by our tentative 
detection of a velocity gradient over the CH3CN source perpendicular to the MMl outflow 
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axis. 

Widespread CO line-wing emission is seen virtually over the entire cluster core area. 
In total, we estimate that a few solar masses of molecular gas have been set into motion. 
This estimate does not include corrections for optical depth effects or outflow gas having 
low radial velocities, however, so the total outflow mass might well be larger by a factor of 
~ 10. The flows are estimated to carry ~16 Mq km s~^ of momentum and a kinetic energy 
of 75 Mq km^s~^, which are lower limits again; correcting the mass for optical depth and for 
CO with Vrad about 0, and correcting the velocities for inclination might well increase these 
values by more than an order of magnitude. 

In clusters of low mass stars, protostellar outflows can impart a significant amount 
of energy and momentum into their surroundings and, perhaps, thereby regulate t he star 



formation rate flNorman fc Silklll980l : iMatzner fc McKed 12000 



Li fc Nakamurall2006[ ). They 



may also play an important role in replenishing turbulence, which hydrodynamic as well 
as magnetohydrodynamic simulations sug gest decays on timescales comparable to a clouds 



free-fall timescale Tg = (37r/32Gp)^" (e.g. lMacLowlll999l ). In order to maintain turbulence 



kinetic energy and momentum has to be provided by the flow at a high enough rate to 
compensate for the loss due to the decay of turbulence. 



The mass of the L1641- N cluster forming core is a few times 100 Mq ( iReipurth et al. 
19981 : iTatematsu e7aDll993h . a certain part of which is in the still inactive 'tail' extending 
south-east of L1641-N/Strom 11; the actual mass involved in directly building the L1641-N 
cluster might be about 100-200 Mp. Us ing a cloud radius of about 0.5 pc (estimated from 
the map shown by lReipurth et al.l (119981 )) and a mass Mcore of 200 Mq yields = 4.1 ■ lO^yr. 
We estimate a one- dimensional velocity dispersion ctid of 0.94 km s~^ for the L1641-N cloud 
core based on the addition (in quadr ature) of the intr i nsic c ore dispersions of 0.72 km s~^ 



and core-core motion, 0.6 km s ^ in ITatematsu et al.l (119931 ): cores 67 and 69). The rate 
Ltmh = Eturh/'Ts at which kinetic energy is lost due to decaying turbulence is then estimated 



to be 0.1 L0 (where 



turb 



3/2M, 



coreC"lD 



The momentum available in the CO flows has the potential to add significant motions 
to the cloud core. If the forward momentum were completely transferred to the 200 Mq 
of core gas, it would accelerate it by at least ~0.1 km s~^ and possibly up to the cloud 
velocity dispersion with the corrections mentioned above. In addition, sideways motions can 
be excited, which does not affect the flows forward oriented momentum. In terms of energy 
supply, even the (lower limit on) L^^ch of 1.21 Lq is greater than Lturb- However, a large 
fraction of that can be expected to be dissipated in shocks. 



The rate Lgain at which the cloud gas gains energy, taking into account radiative losses. 



- 15 - 



can be estimated as follows: we assume that all of the outflow forward momentum is dumped 
in the clump and used to drive the turbulence, and that the outflows accelerate a fraction 
of the clump's mass to a velocity equal to the clump's three-dimensional velocity dispersion 
(TsD- The clump mass set in motion per unit time is then = -Fout/c"3D where F^^t is the 
sum of outflow momentum input rates. Therefore, 



^gam 



1 • Vs 

-Mc(^Ij) = — FoutCTlD 



where ctid is the observed, one- dimensional velocity dispersion. With ctid of 0.94 kms~^ and 
the sum of outflow momentum rates of Font = 13.6 x 10^'^ Mq kms~^ yr^^ from Table 1 the 
energy gain rate is Lgain ~ O.ISLq. This result is comparable to the estimate of the clouds 
rate of energy loss Lturb ~ O.ILq, particularly if keeping in mind that the outflow mass and 
velocity measurements are lower limits due to optical depth and projection, respectively. 
The above calculation indicates that the molecular outflows can indeed provide sufficient 
energy and momentum to sustain cloud turbulence. 

Protostellar outflows are known to be collimated. In regions forming only a small number 
of isolated and well separated protostars, the outflow will only affect a small fraction of the 
total cloud volume and their impact is therefore relatively small. However, in regions of 
clustered star formation, the protostellar density is higher and a larger fraction of the cloud 
volume will be affected. We have found wide-spread CO line-wing emission in L1641-N, 
indicating that the presenc e of several CO outfl ows does affect a signiflcant volume of the 
cluster core. In NGC 1333, lOuillen et al.l (120051 ) found a system of cavities permeating the 
cloud and sugg e sted, if these are relic outflow events, then ~ 10% of the cloud was affected. 



Williams et al.l (l2003l ) found that most of the CO line wing emission in the OMC-2 and 
OMC-3 star forming regions is in broadly dispersed, extended features, which also suggests 
that a large fraction of the cloud volume is affected by outflows. From the small fraction of 
the high velocity emission in young, compact outflows, and using a calculation similar to the 
above, they also estimated that the energy injection rate was sufficient to maintain cloud 
turbulence. 

L1641-N, NGC 1333, and OMC-2/3 are regi ons forming modest sized clusters of a few 
tens to hundreds of stars. In the case of L1641-N, iHodapp fc Deand ( 119931 ) showed that star 
formation has been going on for about 10^ years at 2-3 stars every 10^ years. Assuming 
that protostars drive powerful CO outflows over the flrst 1-2x10^ years of their evolution, 
the 5 CO outflows we found in the region suggest that star formation is continuing at a 
similar pace. NG C 1333 also shows evidence for star formation over million year time scales 
(ILada et al.lll996l). which is still goi ng on, and even the IC 348 cluster might still be in the 
buildup phase (lEisloffel et al.ll2003l ). These modest sized cluster forming regions typically 
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feature a handful to one to two dozens of energetic outflows. Assuming that the formation 
of a more massive cluster (few thousand or ten thousand stars) also proceeds more or less 
steadily over a few million years rather than in one short burst, then one can expect that 
at any time there are many tens or hundreds of stars in their youngest evolutionary stage 
and driving very energetic outflows. It can then be expected that virtually the entire cluster 
volume is affected by outflow activity. 

The action of a population of outflows from low-mass protostars might a lso have conse- 
quen ces for high-mass star formation. In the 'turbulent core' scenario (e.g., iMcKee &: Tan 
20031 ) . massive stars form from high density cores formed within a high pressure, turbulent 
clump. Accretion rates are high enough to overcome the radiation pressure of the central 
object. The formation time scale is comparable to the free fall timescale of the surround- 
ing clump, i.e., to the timescale of turbulence decay. Hence it might not really be crucial 
to maintain the turbulence while the massive star forms. But it might be a prerequisite, to 
keep the turbulence in the high pressure environment until the massive star forming core has 
formed and reached the point when it starts to collapse. 



The 'competitive accretion' scenario (e.g.. lBonnell et al.ll200ll ) requires that gas falls into 
the central region of a forming cluster, allowing the protostars to grow beyond the masses of 
their original envelopes. This infal l might be count e racte d by outflows from the protostars 



in the cluster center. Furthermore, iKrumholz et a 



and collisional growth of protostars (IBonnell et al. 



(120051) claim that competitive accretion 



19981 ) requires efficient removal of kinetic 



energy from the cluster forming gas, which could also be prevented by low-mass protostars 
outflows. 
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Fig. 1. — Overview of the area around the L1641-N cluster. Sohd black contours represent 
1.3 mm dust continuum (IRAM 30 m/MAMBO single-dish). Greyscale: 2.12 fim narrow 
band image, not continuum subtracted. HH objects and SMZ2002 shock excited H2 features 
are marked (e.g., 5-12 is [SMZ2002]5-12). The dashed line marks the boundary of the 
field covered by the wide-field IRAM 30 m/HERA CO map. AN04-10 marks the 1.3 mm 
counte rpart to the AN04 ISO mid-infrared source 10. N31 is source N31 of IChen et al. 
fll993[ ). 
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Fig. 2. — CO emission in the L1641-N cluster area. Blue/red countours: blue-/red-shifted 
CO (IRAM 30 m/HERA). Solid black contours & greyscale: 1.3 mm dust continuum 
(IRAM 30 m/MAMBO single-dish). We suggest that B-E, R-E, and R-SW are part of one 
bipolar outflow driven by the AN04 ISO source 10, with the redshifted lobe being deflected 
by interaction with the L16 41-N giant outflow (R-S). A small bipolar flow (R-SE, B-SE) is 
driven by the N31 source of IChen et al.l (119931 ). R-S and B-N are part of the L1641-N giant 
outflow. B-NE and R-W are resolved into two parallel outflows by our SMA observations 
(see Fig. ED. 
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. 4. — As Fig. [3|, for velocities around the ambient cloud velocity (3.2 to 12.8 km s 
different greyscale stretch). 



5. — As Fig. [3], for redshifted emission (12.8 to 35.6 kms ^). 
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Fig. 6. — 1.3 mm dust continuum emission in L1641-N. Left panel: SMA+MAMBO com- 
bined map (greyscale, blue contours), with the single-dish MAMBO map overlayed as dotted 
contours. The SMA beamsize is indicated as the black ellipse, and the 30% sensitivity ra- 
dius of the SMA primary beam is shown as the dashed circle. Right: the greyscale shows 
a 2.12 /im near- infr ared narrow band image (i^-band continuum + H2 v=l-0 S(l) line, see 



Stanke et al.l (l2002l )). with H2 shocks labeled; dotted contours: ISO 7 /im map (AN04). The 



square marks the position of AN04 ISO source 18. 
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Fig. 7. — C0(2-l) emission from L1641-N. Blue/red contours: blue/redshifted CO; dotted 
black contours: CH3OH. The square marks the position of the AN04 ISO source 18. The 
black/blue elhpses indicate the continuum/CO-hne beam size. Left: greyscale: 2.12/im 
narrow band image (H2 emission + continuum). Black contours: SMA + MAMBO 1.3 mm 
dust continuum (see Fig. [6]). The solid line connects the H2 bow shock Al and the driving 
source MMl. The orientation of the CO outflows clearly deviates from that line. Right: 
low-velocity CO; greyscale, solid black contours: SMA + MAMBO 1.3 mm dust continuum. 
The panels at the bottom show a close-up of the central region. Black lines are drawn to 
emphasize the cavity structures seen at low velocities. 
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Fig. 8. — Distribution of molecular line emission in L1641-N. The greyscale and solid 
black/white contours show line emission integrated from 5 to 8.6 km s~^, the red and blue 
contours show emission integrated from 8.6 to 13.4 km s~^ and from 0.2 to 5 km s~^, re- 
spectively. The dotted contours show the SMA+MAMBO combined continuum map. The 
X symbol marks the ISO 7 /xm peak position. 
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Fig. 9. — Top: CH3CN map (greyscale; integrated over the K = and 1 lines; contours: 
SMA+MAMBO combined continuum map; cross: phase center; eUipse: beam). Bottom: 
CH3CN spectrum of MMl, along with model spectrum (grey line; model parameters: T — 
130K, A^cHgCN = 4 • lO^^cm"^, Av = 4.1kms~"^ source size (y.'42). The dip in the spectrum 
close to the K — 5 line is caused by a number of particularly noisy channels. 
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Fig. 10. — Left: Position of the emission maxima in individual velocity channels for the 
CH3CN (diamonds: K=0; triangles: K=l; squares: K=2; x's: K=3) and OCS lines (+); the 
size of the crosses show the 1-a positional errors as given by the gaussian fitting procedure. 
Offsets are given in arcseconds from the phase center. The square marks the position of the 
MMl continuum peak, the dotted line the orientation of the assumed disk plane (position 
angle 120°). Right: velocity vs. offset from source along the disk plane. The plotting symbols 
are the same as in the left panel. For clarity, velocities in the various transitions are plotted 
with a small offset. Although the positional errors are larger than the positional offsets from 
channel to channel, a trend for increasing velocities (going from southeast to northwest) is 
consistently seen in all transitions. 
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Tablc 1. Parameters of the individual CO outflow lobes. Momentum and kinetic energy 
have been calculated assuming fccn = 6.8 km/s. In the last two rows, mass, -Ekm, and P are 
integrated over the entire area where high velocity CO emission is seen, and includes a 
diffuse component which can not be attributed to individual flows; however, the velocity 
range (close to the ambient cloud velocity) had to be chosen more carefully than for the 
individual CO flow lobes, therefore omitting some of the slow moving line-wing emission, 
-^mech and Fmech in the last two rows are given as the sum over the individual outflow lobes. 
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0.7 


28 
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0.07 
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0.8 
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0.13 
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0.8 
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0.13 
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0.7 
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0.07 


4.4 


0.7 
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0.09 


3.7 
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1.9 
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0.036 


1.4 


0.29 


14.4 


0.04 


2.7 


0.09 


1.1 
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0.044 


0.18 


0.11 
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0.05 
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0.003 


0.1 
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-2.2 . 


..5.0 


0.13 


0.68 


0.38 
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0.18 


24.5 


0.005 


0.16 


R-SW 
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21.8 


0.34 


6.21 


2.0 


12.0 


0.65 


53.0 


0.02 


0.37 


R-S 


9.8 .. . 


26.6 


1.1 


33 


7.6 


14.4 


0.86 
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0.09 


1.3 
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0.95 
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0.68 
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. 39.8 


1.21 


IS 


10.2 








o.m 
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Table 2: Molecular lines detected, and the features where the lines are seen (x means de- 
tected, (x) means marginally detected, and - means not detected). 



I/rest (GHz) 


fobs (GHz) 


Species 


Transition 


MMl 


SW shock 


outflow 


MM3 


219.5604 


219.5601 




2-1 


X 






X 


219.9494 


219.9497 


SO 


6(5)-5(4) 


X 


X 


X 


X 


220.0785 


220.0789 


CH30H 


8(0,8)-7(l,6) E 


X 


X 




(X) 


220.3987 


220.3994 




2-1 


X 




X 


X 


220.6411 


220.6424 


CH3CN 


12(5)-11(5) 


X 








220.7090 


220.7098 


CH3CN 


12(3)-11(3) 


X 






(X) 


220.7303 


220.7307 


CH3CN 


12(2)-11(2) 


X 






(X) 


220.7430 


220.7432 


CH3CN 


12(1)-11(1) 


X 






(X) 


220.7473 


220.7469 


CH3CN 


12(0)-11(0) 


X 






(X) 


229.5891 


229.5875 


CH30H 


15(4,11)-16(3,13) E 


X 








229.7588 


229.7574 


CH30H 


8(-l,8)-7(0,7) E 


X 


X 




X 


230.0270 


230.0261 


CH30H 


3(-2,2)-4(-l,4) E 


X 








230.5380 


230.518 


CO 


2-1 










231.0610 


231.0601 


ocs 


19-18 


X 








231.2208 


231.2185 


i^CS 


5-4 


(X) 








continuum 


(mJy) 






470 






79 



